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Abstract. - 

The clustering of galaxies is well characterized by fractal properties, with the presence of 
an eventual cross-over to homogeneity still a matter of considerable debate. In this letter 
we discuss the cosmological implications of a fractal distribution of matter, with a possible 
cross-over to homogeneity at an undetermined scale Rhomo- Contrary to what is generally 
assumed, we show that, even when Rhomo — > oo, this possibility can be treated consistently 
within the framework of the expanding universe solutions of Friedmann. The fractal is a 
perturbation to an open cosmology in which the leading homogeneous component is the cosmic 
background radiation (CBR). This cosmology, inspired by the observed galaxy distributions, 
provides a simple explanation for the recent data which indicate the absence of deceleration in 
the expansion (q ~ 0). Correspondingly the 'age problem' is also resolved. Further we show 
that the model can be extended back from the curvature dominated arbitrarily deep into the 
radiation dominated era, and we discuss qualitatively the modifications to the physics of the 
anisotropy of the CBR, nucleosynthesis and structure formation. 



One of the most extraordinary findings of the last two decades in observational cosmology 
has been the existence of a network of voids and structures in the distribution of galaxies in 
space. The enormous scales of these structures were completely unsuspected in earlier exten- 
sive observations of galaxy distributions, in which only angular coordinates were measured, 
obscuring the richness subsequently revealed in the third coordinate. These findings have 
become increasingly difficult to reconcile with standard cosmological theories, in which the 
approach to homogeneity at large scales is a central element fm. Observationally, however, not 
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only the scale at which the matter distribution approaches an average density, but the very 
existence of such a scale, remains the subject of intense debate At small scales it is 

well established that the distribution of galaxies is fractal, and the debate can be phrased in 
terms of the deviation from this behaviour towards homogeneity. Some consensus has been 
achieved about the optimal statistical methods to use in the analysis of three dimensional data, 
with disagreement remaining on details of the treatment of some data sets [^[ ^ . The case 
for homogeneity still rests primarily on indirect observational evidence, such as the angular 
data. Independently of the data, however, resistance to the fractal picture is certainly to 
a considerable degree due to the conviction that it is incompatible with the framework of 
the standard theories (see e.g. @), and in particular with the high degree of isotropy of 
the microwave background radiation |8|, |9). In this respect one should note that in standard 
models the origin of radiation and baryonic matter is completely separate, with the latter 
being created in a dynamical process ('baryogenesis') completely distinct from the origin of 
the primordial radiation bath. The isotropy of the latter is therefore not fundamentally tied to 
the distribution of the matter, and the only real constraint is how much any such distribution 
actually perturbs the radiation. In this Letter we do not discuss the details of the evidence for 
or against homogeneity, but rather consider the grounds for these theoretical biases against 
the possible continuation of a fractal distribution to arbitrarily large scales. Our central result 
is that a fractal distribution for matter, even when there is no upper cut-off to homogeneity, 
can in fact be treated in the framework of an expanding universe Friedmann cosmology. 
A fractal JToj, |ll| is a self-similar and intrinsically fluctuating distribution of points at all 



scales, which appears to preclude the description of its gravitational dynamics in the framework 
of the Friedmann- Robertson- Walker (FRW) solutions to general relativity jl| . The problem is 
often stated as being due to the incompatibility of a fractal with the Cosmological Principle, 
where this principle is identified with the requirement that the matter distribution be isotropic 
and homogeneous Q . This identification is in fact very misleading for a non- analytic structure 
like a fractal, in which all points are equivalent statistically, satisfying what has been called a 
Conditional Cosmological Principle [|l0[ |ll|, The obstacle to applying the FRW solutions 
has in fact solely to do with the lack of homogeneity. One of the properties of a fractal of 
dimension D, however, is that the average density of points in a radius r about any occupied 
point decreases as r D ~ 3 , so that asymptotically the mass density goes to zero |h], O. An 
approximation which therefore may describe the large scale dynamics of the universe in the 
case that the matter has such a distribution continuing to all scales is given by neglecting 
the distribution of matter at leading order, relative to the small but homogeneous component 
coming from the cosmic microwave background. We will now show that is indeed a good 
perturbation scheme, and calculate the physical scale characterizing its validity. 

Consider first the standard FRW model with contributions from matter and radiation, for 
which the expansion rate is 

H 2 (t) = = ^( Prad + Pmat ) -k/a 2 (1) 

where a(t) is the scale factor for the expansion, and p la a oc 1/a 4 is the radiation density, and 
Pmat oc 1/a 3 the (homogenous) matter density. The constant k — — H 2 a 2 D (l — VL r — TO ), where 
H a (a G ) is the expansion rate (scale factor) today and Q r (f2 m ) is the ratio of the radiation 
(matter) energy density today to the 'critical' density p c — -^qH 2 . The sign of k determines 
whether the universe is closed (k > 0) or open (k < 0), with k = corresponding to a 'critical' 
spatially flat universe. Given the temperature of the CBR 0, we have^J Q r h 2 « 2.3 x 10 -5 



( 1 )We will neglect here, for simplicity, the minor modifications due to massless or low mass 
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(where h is the Hubble constant in units of lOOMpc/km/s, with a typical measured value of 
h 0.65 JlJ, If we make the simple and natural assumption that galaxies trace the 

mass distribution the value of Sl m depends directly on the determination of the scale of the 
cross-over to homogeneity. If the observed fractal distribution continues to a scale Rhomo, 
above which it turns over to homogeneity, one has 

10 - 3 - D 



tt m = fiio t; (2) 

\ ^homo / 

where fiio is the average density of matter (relative to critical) in a sphere of radius lOMpc/h 
about a galaxy, D is the fractal dimension, and Rhamo is measured in Mpc/h. For Rhomo 
sufficiently large that f2 m < the p m term in (Q) is always sub-dominant, and there is no 
matter dominated era. For simplicity we now consider the limit in which Rhomo - * oo. Solving 
for the scale factor we then have 

a(f) =a (2H o nh)^(l + — ^Hot ) (3) 

which shows how the early time radiation dominated behaviour (a cx fi) changes to the linear 
law a oc t at t w 2iJ~ 1 v / TI7 (red-shift z ~ where 1 + z = a /a). We now discuss 

how in each of these two phases (dominated respectively by the radiation and curvature) the 
fractal can be treated as a perturbation to this solution. When we make numerical evaluations 
below we will use data from galaxy catalogues in Q], which give D w 2 and fiio = 0.007. The 
latter assumes a mass to luminosity ratio of lOh, in solar units, (that estimated for a typical 
spiral galaxy Note that since y/H^ = 0.005h this value requires Rhomo > lOMpc/h for a 

direct transition from radiation to curvature domination. If instead we take the global mass to 
luminosity ratio to be that estimated in clusters (« 300h |l6|]) we require Rhomo > 250 Mpc/h. 

First consider the curvature dominated phase. The radiation is negligible and, at scales well 
within the horizon, we can use Newtonian gravity to describe the solution and its perturbations 
when the self-gravity of the matter is included. The leading solution is simply the free 
expansion of the fractal, with every point moving radially away from its neighbour at a 
constant velocity proportional to its distance i.e. r = H r(t ) = H(t)r(t). To estimate the 
deviation from this flow due to the self-gravity of the fractal, we take a point in the flow and 
integrate the work done against gravity along its trajectory (in the leading order unperturbed 
flow). If the particle moves from an initial position R a , where it feels a total gravitational 
acceleration F(R ), to a final position R — xR , this work done (per unit mass) is simply 
W{R 01 x) — F{R ).R (l — -). The integral is performed along the unperturbed trajectory 
using the fact that the force on the chosen point simply scales as 1/x 2 . A sufficient condition 
for the Hubble flow to apply to a good approximation at all subsequent times is simply that 
W(R ,x — oo) be much less than the kinetic energy of the particle i.e. 

F(R ).R « hilRl (4) 

Noting that the force at the origin of the flow was implicitly taken to be zero, we see that 
the validity of the criterion (|]) will be determined in a fractal by the difference between the 
gravitational force on two occupied points as a function of the distance between them. The 
gravitational force on a point in a fractal has been studied in | p7[ . Its behaviour can be 
understood as the sum of two parts, a local or 'nearest neighbours' piece due to the smallest 



neutrinos, which can easily be incorporated in our analysis. 
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cluster (characterised by the lower cut-off A in the fractal) and a component coming from the 
mass in other clusters. The latter is bounded above by the scalar sum of the forces 

(\F\) < lim [ L ^2^1^r 2 dr ~ L D ~ 2 (5) 
L^ooJ A r l 

so that for D < 2 it is convergent, while for D > 2 it may diverge. If there is a divergence, 
it is due to the presence of angular fluctuations at large scales, described by the three-point 
correlation properties of the fractal. For the difference in the force between two points the 
local contribution will be irrelevant well beyond the scale A, while it is easy to see that the 
'far-away' contribution will now converge as L D ~ 3 , and its being non-zero is a result of the 
absence of perfect spherical symmetry. Noting that as a function of distance R between points 
this component is bounded above by the same behaviour as the force, we write 

(\F(R)\r»a) = A 3 ^± oc R D - 2 (6) 

where the pre-factor A 3 contains non-trivial information about the three-point correlation 
function of the fractal. These convergence properties of the relative force on two points are 
enough to draw a simple conclusion from the criterion (jij): For a fractal in Hubble flow there 
is always a scale above which its evolution will be well described by continued Hubble flow for 
all subsequent times. 

We now apply this to the Universe, and estimate the physical scale today R a up to which the 
unperturbed 'no matter' Hubble flow can be maintained right through the curvature dominated 
era. Given that this era begins at a redshift z « l/yO^ ~ 200/i, we require 

F(R)R\ z =200h w {200h)F{R o )R o < ^H 2 R 2 a -> R > 20n 10 (200/i)A 3 (7) 

What is observed is Hubble flow with deviations (peculiar velocities) only at 'cluster' scales 
~Mpc. Taking our estimate for SIiq we thus require A 3 < \/2i)h. A fractal with a very 
weak three point correlation is one which has a very isotropic angular projection, so if the 
Universe is indeed a fractal a small value of A3 would be expected. It is simple also to derive 
an expression for the peculiar velocity (small compared to the Hubble flow velocity Vh) which 
implies a simple linear relation, just as in standard perturbed homogeneous cosmology [ [b8[ , 
between the local force and the velocity perturbation (Av/vh)(R) oc (F(R)/R). This relation, 
for which there is apparently observational support, is usually used to determine an unknown 
constant (the 'bias' factor) Jl8| |. In the present framework it can in principal be used to extract 
information about the total mass density and the constant A3, which in turn can be related 
to angular data (and ultimately measured directly in forth-coming red-shift surveys). 

Here we have assumed that the fractal extends to arbitrarily large scales. For a finite Rhomo 
the analysis can be easily modified, by breaking the integrals at the appropriate scale. The 
part from scales greater than Rhomo will give a contribution which can be reabsorbed within 
the Hubble flow, while the perturbation will maintain the same scaling at smaller scales. The 
central point which we emphasize is that the scale R is only indirectly related in this case to 
the homogeneity scale, and remains finite as Rhomo - > 00. 

We have thus seen that an open FRW universe is always a good approximation beyond 
some finite scale if matter is distributed as a simple fractal up to an arbitrarily large scale. In 
particular such an open model - because it is dominated by the kinetic energy of the Hubble 
flow - can explain naturally how large structures can co-exist with almost perfect Hubble 
flow. We further note a few other of its striking features: (i) Since the Universe is to a good 
approximation in completely free expansion at large scales with a(t) oc t, we have deceleration 
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parameter q Q sa 0. This is a good fit to recent supernovae observations |2C|| . Rather than being 
due to the effect of an unknown 'anti-gravitational' component which mysteriously cancels 
the decelerating effect of the matter on the expansion the effect is due to the decay 
towards zero of the matter density on such scales, (ii) The expansion age of the Universe is 
t = H^ 1 w lO/i -1 ss 15 billion years, larger by 50% than in the standard matter dominated 
case. This value is comfortably consistent with the estimated age of globular clusters (the 
oldest known astrophysical objects) f 1.5 ± f .3 billion years pl|] . (Hi) The size of the horizon 
today is -Rff(i ) ~ — ^cH^ 1 \nfl r ss 20, OOOMpc/h, a factor of about three larger than in the 
standard case. 

So far we have considered the model only in the curvature dominated era i.e. back to 
red-shift z m 200h. For this last point above however we have extrapolated the model back 
to the radiation dominated era, assuming that the effect of the matter distribution can also 
be consistently treated as a correction in this epoch to the FRW solution without matter. We 
now justify this assumption, and then discuss some of its consequences for a specific cosmology 
of this type. For the former we simply treat the fractal as a set of perturbations to the energy 
density in a manner analogous to the way such perturbations may be treated in the standard 
framework. There the criterion one would use to apply the uniform Hubble flow to describe 
the growth of the horizon is simply that such perturbations be small at the horizon scale i.e. 
^j\hor be small (where p is the homogeneous energy density i.e. that in the radiation). In a 
fractal the perturbations are non-analytic and Sp has no meaning as defined in the standard 
case. We can however write down the mean mass (or energy) at the scale of horizon about an 
occupied point. Taking this as the appropriate Sphor is clearly the right adapted criterion, as 
the fractal is simply made of voids and structures, and voids clearly will not perturb the flow. 
We thus require that the fractal obey 



where for p ma t we have taken the mass inside the (comoving) horizon Rh(z) at red-shift z. In 



Rh ~ 600 Mpc so that, for a fractal with D = 2, we have at this red-shift 8h ~ ^io- In order 
that the fractal matter be indeed a small correction at this red-shift we require, approximately, 
Oio < f , which holds comfortably even if there is much more dark matter than we assumed in 
obtaining our estimate f^io ~ 0.01 . We can thus continue to use the FRW solution back to an 
arbitrary red-shift for a fractal distribution of matter extending to the corresponding scales, 
provided that the condition (||) holds. 

To make a link with central observations in cosmology such as the microwave background 
and nucleosynthesis, we need to specify a precise model. In the spirit of this letter we now 
consider here a radical (but very simple) possibility for a cosmology which makes use of the 
results we have presented: We consider a universe which at very early times (deep in the 
radiation dominated era) is a radiation bath at a given temperature with superimposed fractal 
perturbations in baryon number up to the arbitrarily large scale Rhomo, and down to some 
scale A (see below). Note that positing a very different distribution for matter and radiation 
does not represent a loss of simplicity in comparison to standard models, which generically 
envisage the (almost) homogeneously distributed matter as coming from a dynamical process 
('baryogenesis') completely distinct from the origin of the primordial radiation bath. Instead 
of fixing the initial condition on a homogeneous baryon to photon ratio (ng/ti 7 ~ 10~ 9 ), with 
some independent superimposed spectrum of analytical fluctuations, we specify our fractal in 
baryon number from the properties observed in the distribution of visible matter at large scales 
today. This is the appropriate normalization given that these perturbations are simply frozen 




(8) 



the radiation dominated era (for z > l/yTV) we have Rh(z) 




£=. At z — 10 3 this gives 
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at all but very small scales in the curvature dominated era. In particular we take D — 2 and 
the normalization of the mass given by £lio- Interestingly, stated in terms of the parameter 5h 
above, these values correspond to the special case that 8h is constant, and of order one. Below 
the lower cut-off scale A we take the distribution to be smoothed, with the corresponding 
density f2io(10/A) (where A is the comoving scale given in units of Mpc/h). A natural lower 
bound for this scale is that characterizing the baryon diffusion (up to the corresponding time), 
which will smooth any inhomogeneity at smaller scales. 

The fluctuations in the temperature of the microwave background depend on the 'intrinsic' 
fluctuations imprinted at time of decoupling of the photons, plus the fluctuations induced 
in their propagation from their last scattering. We consider here only the former as they are 
typically the dominant effect for perturbations which are essentially frozen on the scales we are 
interested in. Relative to the standard critical or near critical mass universe there are two main 
features to note here. First, photon decoupling will be modified greatly: While in the standard 
case there is a global baryon density which determines the time/temperature of decoupling, 
here the relevant baryon density varies enormously - for a photon in a void it is zero, while for 
one in a structure it is the local density of baryons associated with the lower cut-off scale A, 
a density which can be several orders of magnitude greater than in the standard case (if we 
take A as small as the baryon diffusion distance). However, since the decoupling temperature 
is only logarithmically sensitive to this parameter, the decoupling of photons in structures will 
still occur around redshift z ~ 10 3 . On the other hand, if the scale Rhomo is so large that 
there are voids of the order of the horizon scale at this time (at z = 10 3 we found Rh ~ 600 
Mpc), most photons will decouple at the much earlier time of electron-positron annihilation 
since after this time they find themselves in a neutral environment. Thus the 'thickness' of 
the surface of last scattering will be very much greater than in the standard case, essentially 
consisting of two stages of 'void decoupling' (at a redshift of ~ 10 9 ) and 'structure decoupling' 
at redshifts comparable to the standard one. 

The other main difference relative to the standard case is that, because of the extremely 
low background density in the model, the effect of the hyperbolic geometry is much greater. In 
particular, at high red-shift (z^/TT^ >> I), the angle 9 subtended by a scale of given physical 
size Xg is 9 « Ag/ H~ 1 v / H7 which means that the physical scale corresponding to 1° on the sky is 
I0 4 Mpc i.e. of the order of the horizon scale today, and a factor of about 100 larger than in the 
standard case. Considering the possibility that the fractal extends to such enormous scales, we 
can make a naive estimate of the amplitude of fluctuations on the microwave sky: Adapting the 
Sachs- Wolfe formula as derived for the standard case of analytical fluctuations, the maximum 
amplitude of temperature variation between two photons subtending some angle should be 
estimated by the energy density fluctuation represented by a structure at the corresponding 
physical scale. At 10° on the sky we then have 



where Zd is the decoupling red-shift. Thus, the largest effect will come from the photons which 
decouple last, for which this maximum amplitude will be ~ 10~ 5 , which is comparable to 
the average amplitude observed at this scale by COBE ||. On the other hand, for a modest 
value of the crossover scale to homogeneity (e.g. Rhomo ~ lOOMpc/h), the effect of the 
fractal distribution of matter on the photons at decoupling will only be visible at scales much 
smaller even than those which will be probed by the Planck satellite mission p2|. A detailed 
study of the perturbations induced by propagation of photons through an expanding structure 
of this kind - requiring techniques quite different to the standard ones treating analytical 
perturbations - will be required to see if it is possible to produce perturbations at the levels 




1 Pb(Aiq) = I fl B (X ) 

3 Pradn(zd) 1 + Zd 3£l r 



(9) 
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observed by experiment at the angles which have been probed to date. 

Finally a few brief comments on nucleosynthesis and structure formation. The cooling rate 
of the plasma is the same as in the standard case, and the results of nucleosynthesis will 
depend on the local baryon to entropy ratio, which as discussed above is related to the scale 
A. If this scale is larger or comparable to the horizon scale at nucleosynthesis the amount of 
helium produced will not differ much from the standard case - it is essentially independent of 
the baryon to entropy ratio - while the residual densities of deuterium etc. will be lower. For 
a smaller (i.e. sub-horizon) A the effect of inhomogeneities will be important, but a reliable 
calculation for the effect becomes very difficult to perform. It is unlikely however to modify 
the tendency for lower values of the 'trace' elements, since this arises due to the fact that 
the elements are synthesized in denser regions compared to the standard case. If, on the 
other hand, the fractal (up to a finite sub-horizon Rhomo) is formed after nucleosynthesis, the 
appropriate value for the local density would be that at Rhomo- We note that for modest values 
of this scale {Rhomo ~ 50 — WOMpc/h) this would correspond to standard nucleosynthesis if 
there is ratio of dark to visible baryons comparable to the mass to light ratio inferred in 
clusters Clearly some physics quite different to that at work in standard models would 
be required to make it possible to generate such a structure between nucleosynthesis and the 
curvature dominated era when (as we have noted) the structure gets frozen in at all but small 
scales. These and other considerations - in particular a more detailed study of the microwave 
background fluctuations in these kinds of structures - we will pursue in forthcoming work. 
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